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Abstract 

The spectrum of quarter BPS dyons in M = 4 supersymmetric string theories can 
change as the asymptotic moduli cross walls of marginal stability on which the dyon can 
break apart into a pair of half BPS states. In this paper we classify these marginal stability 
walls and examine this phenomenon in the context of exact dyon spectrum found in a class 
of M = 4 supersymmetric string theories. We argue that the dyon partition functions 
in different domains separated by marginal stability walls are the same, but the choice 
of integration contour needed for extracting the degeneracies from the partition function 
differ in these different regions. We also find that in the limit of large charges the change in 
the degeneracy is exponentially suppressed compared to the leading contribution. This is 
consistent with the fact that in the computation of black hole entropy we do not encounter 
any change as the asymptotic moduli fields move across the walls of marginal stability. 
Finally we carry out some tests of S-duality invariance in the theory. 
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1 Introduction and Summary 

In a series of papers we computed the exact degeneracy of quarter BPS dyons in a 
class of M = 4 supersymmetric string compactifications in an appropriate corner of the 
moduli space of these theories [1-3] verifying and generalizing earlier conjectures [4-6]. 
Alternative approaches to this problem leading to similar results have also been devel- 
oped [4,7-10]. The result for the degeneracy takes the form of integration over a three real 
dimensional subspace (a contour) of the Siegel upper half plane parametrizing genus two 
Riemann surfaces, and the integrand involves inverse of a certain meromorphic modular 
form of a subgroup of the Siegel modular group. 

It is well known however that for a dyon with a given set of charges the moduli space 
of M = 4 supersymmetric string theory contains subspaces on which the original dyon 
becomes marginally unstable against decay into a pair of other dyons [11,12]. In particular 
the moduli space contains walls of marginal stability - codimension one subspaces - on 
which the mass of the quarter BPS dyon becomes equal to the sum of masses of a pair of 
half BPS dyons whose charges add up to that of the original quarter BPS dyon. On this 
subspace the original dyon becomes marginally unstable against decay into this pair of 
half BPS dyons, and typically the spectrum of the original dyon changes discontinuously 
as we move through these marginal stability walls in the moduli space [11,12]. Thus 
an important question is: how does the dyon spectrum computed in [1-3] change as we 
move away from the particular corner of the moduli space in which the degeneracy was 
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computed? 

A glimpse of this issue was already seen in the analysis of [1-3] where it was found 
that even in the corner of the moduli space where the result was computed, - in a weakly 
coupled type IIB string theory compactified on a certain orbifold, - the result for the 
degeneracy changes discontinuously as the angle between certain pair of circles of the 
compact manifold passes through zero. This change could be attributed to the fact that 
precisely at this point the system under consideration became marginally stable against 
decay into a pair of half BPS states. However the change was such that the expression 
for the degeneracy continues to be given by a similar integral with identical integrand, 
but the contour over which the integral is to be performed gets changed. If we try to 
deform the new contour into the original contour we encounter a pole of the integrand 
and hence the two contributions differ by the residue at the pole. One can also interpret 
the result by saying that the dyon partition function formally remains the same as we 
move through the marginal stability wall, but the point in the Siegel moduli space around 
which we should series expand the partition function to extract the degeneracies changes 
as we move through these walls. 

In this paper we classify these marginal stability walls in the moduli space of the 
theory and explore what happens when we move across these walls. Although the walls 
are complicated codimension one surfaces in the moduli space, one gets a simpler picture 
by regarding them as curves in the axion-dilaton moduli space for fixed values of the other 
moduli. We find that these curves are circles and straight lines, and could intersect on 
the real axis or at zoo, but have no intersection in the interior of the upper half plane. 
Furthermore although the slopes of the lines and the radii and the centres of the circles 
depend on the charges and other moduli, the points where they intersect the real axis 
and each other are universal. As a result a region bounded by these curves has universal 
vertices but boundaries which depend on the other moduli and charges. One such region 
has been displayed in Fig{T] in £j3l As we move from one of these regions to another, we 
cross marginal stability walls and as a result the dyon spectrum could change. 

Typically most of these walls lie outside the domain in which the approximation made 
in our computation of the degeneracy can be trusted. However using T- and S-duality 
invariance of the theory we can extract useful information about these walls. A useful 
input in reaching this conclusion is the observation made in [13] that under an S-duality 
transformation the degeneracy formula does not remain invariant, but there is a change in 
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the integration contour. This can be attributed to the fact that a duality transformation 
acts both on the charges as well as the asymptotic values of the moduli fields, whereas 
the degeneracy formula is computed for different charges but in the same region of the 
moduli space, in a sense that will be made precise in §H The apparent lack of duality 
symmetry of the degeneracy formula is due to the fact that the region in which the 
degeneracy formula is calculated and the region obtained after a duality transformation 
are typically separated by walls of marginal stability. Thus the knowledge of how the 
contour changes under a duality transformation can be used to extract information about 
how it changes from one region of the moduli space to another as we move through these 
marginal stability walls. This way the change in the spectrum across any wall of marginal 
stability can be encoded as the result of changing the integration contour leaving the 
integrand unchanged, or equivalently as shifting the point around which we expand the 
partition function to extract the degeneracy. 

One can try to find a physical interpretation of these changes in the degeneracy by 
explicitly evaluating the residues at the poles picked up by the contour as we move it from 
the position associated with one side of a wall of marginal stability to that on the other 
side of the wall. We can do this easily in the context of the marginal stability wall found 
in the original analysis of [1]. We find that the change in the degeneracy is proportional to 
the total degeneracy of the half-BPS electric and magnetic states into which the original 
dyon can decay on this particular marginal stability wall. 

Given these results, it is natural to ask what happens in the large charge limit where 
the statistical entropy given by the logarithm of the degeneracy can be compared with 
the black hole entropy. We find that in this limit the change in the degeneracy as we 
move across a wall of marginal stability, - encoded in the residues at the poles which 
we encounter while deforming the new contour to the old one, - is exponentially sup- 
pressed relative to the leading contribution to the degeneracy. As a result the change in 
the statistical entropy as we move across the wall of marginal stability is exponentially 
suppressed. This is consistent with the fact that the entropy of the corresponding black 
hole is controlled solely by the charges and is independent of the asymptotic moduli due 
to attractor mechanism. Thus as the asymptotic values of the moduli fields move across 
the wall of marginal stability the answer for the black hole entropy does not change. Our 
results indicate however that if we are able to incorporate non-perturbative (in inverse 
charges) corrections in the computation of the black hole entropy then we should see a 
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dependence of the entropy on the asymptotic values of the moduli fields, possibly along 
the line of [14,15]. 

Since we use duality invariance to find how the contour should be deformed as we 
move across various walls of marginal stability, one could ask if there is any non-trivial 
test of duality that one could perform. If we could identify duality transformations which 
leave a region invariant, - in a sense that will be made precise in §HJ - then the expression 
for the degeneracy should not change under such a duality transformation. This requires 
that under such a duality transformation either the contour should remain unchanged or 
it should move to a new position such that in deforming it from the new to the old position 
we do not encounter a pole in the integrand. This can then be explicitly tested. There 
is also a possibility that a duality transformation maps one region to another such that 
the degeneracy in each of these regions can be computed directly. In this case we have an 
independent result of how the contour should change as we move from the first region to 
the second and this can then be compared with the predictions coming from duality. We 
identify some duality transformations of these types and carry out the required consistency 
checks. 

The rest of the paper is organized as follows. In §|2] we review the results of [1-3] 
about the degeneracy of quarter BPS dyons in a class of M = 4 supersymmetric string 
theories. In §|3] we determine the locations of the walls of marginal stability in TV = 4 
supersymmetric string theories. In §H] we determine how the different domains of the 
moduli space, bounded by the marginal stability walls, are mapped to each other under 
T- and S-duality symmetries of the theory and use this information to determine how 
the spectrum should change as we pass through a particular marginal stability wall. In 
§[5] we show that the change in the statistical entropy as we move across the marginal 
stability walls is non-leading compared to the full entropy. In §|6]we perform some tests 
of S-duality invariance of the theory. 

Finally we would like to remind the reader that our analysis will focus on the marginal 
stability walls associated with decay of the dyon into a pair of half-BPS states. It will be 
worth exploring if there are interesting phenomena associated with decay into a pair of 
quarter BPS states. 

Some related issues have been discussed in [13]. 
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2 Dyon Spectrum in a Class of J\f = 4 Super symmet- 
ric Models 



In this section we shall review the dyon spectrum in a class of Af = 4 supersymmetric 
models analyzed in [1-3,5,6]. Somewhat different approaches leading to similar results 
have been developed in [4,7-10]. 

These theories are constructed by taking a orbifold of type IIB string theory on 
MxS 1 xS 1 where M is either K3 or T 4 . The generator g of the TL^ group involves 1/N 
unit of shift along the circle S 1 together with an order N transformation g in M. g is 
chosen so that it commutes with an Af = 4 supersymmetry algebra of the parent theory. 
Thus the final theory has Af = 4 supersymmetry. 

The description of the theory given above will be referred to as the first description of 
the theory. Another useful description is obtained by a series of duality transformations. 
We first make an S-duality transformation in the type IIB theory. Next we make an 
R — > 1/R duality on the circle S 1 that takes type IIB string theory on M x S 1 x 5 1 
to type IIA string M x S 1 x S 1 where S 1 is the circle dual to S 1 . Finally using six 
dimensional string-string duality we relate this to a heterotic string theory on T 4 x S 1 x S* 1 
for M = K3 and type IIA string theory on T 4 x S 1 x S 1 for M = T 4 . Under this duality 
the transformation g gets mapped to a transformation g that acts only as a shift on the 
right-moving degrees of freedom on the world-sheet and as a shift plus rotation on the 
left-moving degrees of freedom. In the final theory, obtained by taking the orbifold of 
heterotic or type IIA string theory on T 4 x S* 1 x S* 1 by a 1/N unit of shift along S* 1 
together with the transformation g, all the space-time supersymmetries come from the 
right-moving sector of the world-sheet. We shall call this the second description of the 
theory. 

These theories typically have several moduli fields. In the second description one 
complex modulus scalar arises from the axion (a) - dilaton (0) combination, where the 
axion by definition is the scalar field obtained by dualizing the NSNS sector 2-form field. 
We shall denote by r the combination a + iS where S = e~ 2li> . Other moduli fields arise 
from the g invariant components of the metric, antisymmetric tensor field and gauge fields 
(in case of heterotic string theory) along the six dimensional internal torus and may be 
encoded in an r x r matrix valued field M satisfying 

MLM T M. M T = M , (2.1) 
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where L is a matrix with 6 eigenvalues +1 and (r — 6) eigenvalues —1. Here r is the 
rank of the gauge group and depends on the specific model being considered. These U(l) 
gauge fields arise from the g invariant ten dimensional gauge fields (in case of heterotic 
string theory) as well as g invariant components of the metric and the 2-form field with 
one leg along the internal torus and one leg along the non-compact directions. 

Following the chain of dualities relating the first description to the second description 
one can work out the origin of the various fields in the first description. In particular the 
modulus r = a + iS can be shown to correspond to the complex structure modulus of the 
torus spanned by the S 1 and S 1 directions. 

We shall now give the precise relationship between some of the moduli fields and the 
geometric quantities associated with the compactification in the second description and 
also give precise expression for some of the charges in terms of physical quantum numbers 
carried by a state. To do this we need to fix a normalization convention for the coordinates 
along S* 1 and S 1 . We choose x 4 and x 5 to be the coordinates along S* 1 and S* 1 respectively 
and choose their periods before orbifolding to be 2-Ky/a' and 2nN\fa' respectively. Thus 
after orbifolding both can be regarded as having period 2it\fa'] however since periodicity 
under a 2n\^a' translation action S* 1 also involves an action of the generator g of the 
internal symmetry group, the momentum along S 1 , measured in units of l/x/a 7 , 
could be fractional, - some multiple of 1/N, - if g acts non-trivially on the state. A 
string will be said to carry one unit of winding along S* 1 (or S* 1 ) if, as we go once around 
the string, the x 4 (x 5 ) coordinate shifts by 2-Ky/a' '. Thus an untwisted sector state whose 
coordinate along S* 1 changes by multiples of 27riV\/a 7 will carry winding charge along S 1 
in multiples of N, but twisted sector states can carry generic integer winding charges. A 
single H-monopole associated with S 1 will correspond to an array of NS 5-branes wrapped 
on S* 1 x T 4 and placed at intervals of 27ry / a 7 along S 1 . Finally the original Kaluza-Klein 
monopole associated with S 1 , represented by a Taub-NUT space with an asymptotic circle 
of radius 2nN\fa i along x 4 , will develop a "M N singularity at its centre after the orbifolding 
and has to be regarded as carrying N units of Kaluza-Klein monopole charge associated 
with S 1 . Thus the Kaluza-Klein monopole charge associated with S* 1 will be quantized in 
units of N. Similar definition can be given for the H and Kaluza-Klein monopole charges 
associated with S 1 , but in this case the normalization is straightforward and both the 
charges are allowed to take arbitrary integer values. 

Let x^ (0 < ji < 3) denote the coordinates along the non-compact coordinates. For 
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our analysis it will be useful to study in detail a subsector of the theory in which we 
include only those gauge fields which are associated with the 4/i and 5/x components of 
the metric and the anti-symmetric tensor field, only those components of M which encode 
information about the components of the metric and the anti-symmetric tensor field along 
S 1 x S 1 , the axion-dilaton field, and the four dimensional metric. In this subsector there 



are four gauge fields AjL (1 < i < 4) and a 4 x 4 matrix valued field M satisfying 



M 1 = M. MLM 



L = 




h 



h 




(2.2) 



The fields Af and M are related to the ten dimensional string metric Gmn and 2-form 
field B MN via the relations [16,17]: 



G m 
M 



mn i 
t-1 



= <7(io) 
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G 

BG- 1 
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J ->mn — -"mn ) 

G- l B 
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m, n 



4,5 



-XG- 



A(m-1) 



IR(IO) _ R A{m-3) 



4<m,n<5, < fi,u < 3 . 



(2.3) 



A general dyonic state in the theory is characterized by an r dimensional electric 
charge vector Q and an r dimensional magnetic charge vector P. However if we consider 
a dyon that is charged only under the gauge fields (1 < i < 4) introduced in (12. 3p . 
the corresponding charge vectors can be taken to be four dimensional. If we consider a 
state with momentum n and winding — w along S 1 , momentum n' and winding — w' along 
S 1 , Kaluza- Klein monopole charge and H-monopole charge — W associated with S* 1 
and Kaluza-Klein monopole charge N' and H-monopole charge —W associated with S l , 
then we define the four dimensional electric charge vector Q and magnetic charge vector 
P characterizing the state to be 



Q 



n' 
w 
\w'J 



p 



(W\ 
W 
N 
\N'J 



(2.4) 



From our earlier discussion it follows that n' is quantized in units oil/N, N' is quantized 
in units of N and all other quantum numbers appearing in (12.41) are quantized in integer 
units. The precise relation between Q, P and the electric and magnetic charges associated 
with the gauge fields A$ has been derived in [18] by working in the a' = 16 unit. 
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Restricted to this subspace the T-duality transformation of the theory is parametrized 
by a 4 x 4 matrix Q preserving the charge lattice and satisfying 

n T m = L. (2.5) 

The action of T-duality takes us from a charge vector (Q,P) to another charge vector 
(Q',P') given by 

P'=(Q r )- 1 P, Q / = (Q T )- 1 Q, (2.6) 
and the moduli field M to M' given by 

M' = VLMn T , (2.7) 

leaving r unchanged. From eqs. fl2.5p . (12.61) it follows that the following inner products 
are invariant under T-duality transformation 

Q 2 = Q T LQ = 2(nw + nV), P 2 = P T LP = 2(NW + N'W), 

P-Q = P T LQ = Nn + Ww + N'n' + W'w' . (2.8) 

These combinations are independent of the moduli fields M. Using the moduli field we 
can construct more general T-duality invariant combinations like Q T MQ, P T MP and 
Q 7 'MP. We shall make use of these quantities in £j3j 

In the full theory Q and P are r dimensional vectors and Q and M are r x r matrices, 
and L is also an r x r matrix with 6 eigenvalues 1 and (r — 6) eigenvalues —1. However 
the form of eqs.f l2.5p -f l2.7l) remains the same. 

The theories under consideration also have S-duality symmetry which leaves the field 
M unchanged, changes the vector (Q, P) to another vector (Q", P") via the formula 

(?)"("?)(?)• i°s)^- ™ 

and transforms r = a + iS to 

r" = ^, (2.10) 

where the Ti(iV) group of matrices is defined by the conditions 

a5-pj = l, a,5 = I mod N, 7 = mod N, (3 E TL . (2.11) 

We now consider in the first description of this theory a configuration with a single 
D5-brane wrapped on Ai x S 1 , Q\ Dl-branes wrapped on S 1 , a single Kaluza-Klein 
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monopole associated with the circle S 1 , momentum —n/N along S 1 and momentum J 
along S 1 [7]. Since a D5-brane wrapped on M. carries, besides the D5-brane charge, — ( 
units of Dl-brane charge with ( given by the Euler character of M. divided by 24 [19], the 
net Dl-brane charge carried by the system is Q\ — Q. By following the chain of dualities 
described earlier and a suitable sign convention for the charges in the first description, 
we can map this to a configuration in the second description with momentum —n/N 
along S 1 , a single Kaluza-Klein monopole associated with S 1 , {—Qi + C) NS 5-brane 
charge wrapped along T 4 x S 1 , J NS 5-brane charge wrapped along T 4 x S 1 and unit 
fundamental string winding charge along S l [1]. In particular the Kaluza-Klein monopole 
charge associated with S 1 in the first description gets mapped to the fundamental string 
winding number along S 1 in the second description and the D5-brane wrapped on Ai x S 1 
in the first description gets mapped to Kaluza-Klein monopole charge associated with S 1 
in the second description. Using (12. 4p we see that this corresponds to the charge vector^] 



Q 



( o \ 

-n/N 


-1 



(Qi-C\ 
-J 
l 
o 



(2.12) 



This gives 



Q 2 = 2n/N, P 2 = 2(Q 1 -(), Q ■ P = J . 



(2.13) 



The counting of states of this system in the weak coupling region of the first description 
was carried out in [1-3]. We shall now summarize the results of this analysis. We denote 
by d(Q, P) the number of bosonic minus fermionic quarter BPS supermultiplets carrying 
a given set of charges (Q, P), a supermultiplet being considered bosonic (fermionic) if it is 
obtained by tensoring the basic 64 dimensional quarter BPS supermultiplet, with helicity 
ranging from — | to |, with a supersymmetry singlet bosonic (fermionic) state. Our result 
for d(Q, P) is 



d(Q,P) = - 



dp da dve 



iri{NpQ 1 +aP A /N+Z5Q-P) 



(2.14) 



where C is a three real dimensional subspace of the three complex dimensional space 



1 Recall that an NS 5-brane wrapped on T A x S* 1 represents an H- monopole associated with S 1 and an 
NS 5-brane wrapped on T 4 x S 1 represents an H-monopole associated with S 1 . 
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labelled by (p = p% + ip 2 , a = ai + ia 2 , v = V\ + iv 2 ). C corresponds to the subspace 



p 2 = Mi, a 2 = M 2 , v 2 = -M 3 , 
< pi < 1, < ai < N, < vi < 1 , (2.15) 

Mi, M2 and M3 being large but fixed positive numbers with M3 << Mi, M2, and 

$(p,a,w) = e 2 ™(«p+^+v) 

1 AT-l , ^ 

|| || (1 - exp (2m{k a + Ip + jw))) s=0 6 . 

fc',;>o,j<oforfc'=;=o 

(2.16) 

The coefficients c^' 8 \u), 5, 7 encode information about the spectrum of two dimensional 
superconformal u-model with target space Ai and are defined as follows. First we define 

F {r ' s \r, Z) = ^Tr RR -gr ^(_l)^+^ e 2^rLo e -2^Lo e 2«F L ^ ? 0<r, s<N — 1, 

(2.17) 

where Tr RR .gr denotes trace over all the Ramond-Ramond (RR) sector states twisted 
by g r in the SCFT described above, L n , L n denote the left- and right-moving Virasoro 
generators and Fl and Fr denote the world-sheet fermion numbers associated with left 
and right-moving sectors in this SCFT. In defining L and L of a state we subtract 
cx,/24 and cr/24 from the conformal weights of the corresponding operators so that the 
RR sector ground state has Lq = Lq = 0. Due to the insertion of (— 1) Fr factor in the 
trace the contribution to F^ r ' s ^ comes only from the Lq = states. As a result F^'^ does 
not depend on f. Furthermore, using the existence of an SU{2)i J x SU(2)r R-symmetry 
current algebra in this theory one can show that F^ r ' s \r, z) have expansions of the form 

1 

F (r,s)( TjZ) = J2 4 r ' s) (4n-j 2 )e 2 ™ r+2 ^ 2 , (2.18) 

6=0 je2'Hj + b,n£'Hj/N 
4n-j 2 >-b 2 

2 Thc sign of V2 chosen here differs from the ones used in [1-3]. As will be discussed later, this choice 
of sign is valid for a specific choice of sign of the axion field; for the other choice the sign of V2 needs to be 
reversed. The other difference is that we have put the M3 << Mi, Mi condition to make explicit the fact 
that while interpreting d(Q, P) as the coefficients of the Fourier expansion of l/$(/0, ct, v) as in eqs. (|2.2f p . 
(I2.22p we need to first carry out the expansion in powers of e 2lTl P and e 27Tl,J / N and then carry out the 
expansion in powers of e~ 27 ™. This was implicit in the results of [1-3] but was not stated explicitly. 
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for some coefficients c£ («). This defines the coefficients c^'"'(u). We now define 



Q r , s = N ^)(0) + 2cf' s) (-f)) 



N-l 



a: 



24N 



Qo '° ~ 2N ^ Q °' s (T 



-2iris/N 



-2nis 



s=l 



/Ny 



7 



24N 



Qo,o ■ 



(2.19) 
(2.20) 



This defines all the coefficients appearing in (12.141) . 

As an alternative to ( 12.141) . (12.151) we can express d(Q,P) as 

d(Q,P)=g(^Q 2 ,^P 2 ,Q-P 



(2.21) 



where g(m,n,p) are the coefficients of Fourier expansion of the function l/$(p, a,v): 

1 



g(m,n,p) e 



2iri(m p+n a+pv) 



Let us denote by G the group of 4 x 4 matrices generated by: 



(2.22) 



A B 
C D 



A B 
C D 



A B 
C D 



(a 
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o\ 







1 










c 





d 
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\o 








1/ 




/ 


1 







o\ 


-1 
























1 









-1 


0/ 


(1 








/' 


\ 




A 


1 


/' 















1 


-A 












1 


/ 





ad — be 



c = mod N, a, d = 1 mod AT 



A, e ^, 



(2.23) 



and let G denote the group of 4 x 4 matrices satisfying the requirement that: 

g e G iff U- l gU eG, 

where ^ 



U 



( 

-Va 7 o 

\[N 
\ 







-l/V^V 



-l/ v / iV\ 



/ 



(2.24) 



(2.25) 



3 Even though U contains factors of y/N the elements of G are integers due to the fact that for 
A B 
C D 



G G, the elements of C are multiples of N [5]. 
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A. B \ ~ 

An element ( ~ ~ J G G induces a transformation on (p, a, v) via the relations: 

I I) = (An + S){dn + by\ n=(£ fj , (2.26) 

and one can show that $ transforms cLS cL Sic gel modular form of weight k under this 
transformation [3]: 

a, 5*) = det{CQ + D) k $(p, a, u) , (2.27) 

where 

iV-l 

fc = ^E c o M (°)' ( 2 - 28 ) 

This finishes our review of the main results. However two points about the degeneracy 
formula given above need special mention. Eqs. fl2.14p and (I2.2ip are equivalent only if 
the sum over m, n, p in (I2.22p are convergent on the contour C. For the choice of C given 
in ( I2.15P this requires that $(p, cr,^) -1 has a power series expansion in positive powers 
of e 2m ? and e 2 " Kla l N except possibly for a finite number of negative powers, and that the 
coefficient of any given term in this double power series expansion has an expansion in 
positive powers of e -2 ™ except possibly for a finite number of negative powers. This in 
particular requires that the sum over m and n in (12. 22ft are bounded from below, and 
that for fixed m and n the sum over p is bounded from above. By examining the formula 
(12.161) for $ and the result that the coefficients c£ (u) are non-zero only for u > — b 2 , we 
can verify that with the exception of the contribution from the k' = I = term in this 
product, the other terms, when expanded in a power series expansion in e 2mp , e 2mu and 
e 2mv , does have the form of (12.221) with p bounded from above (and below) for fixed m, n. 
However for the k' = I = term, which gives a contribution e~~ 2mv / (1 — e~ 2mv ) 2 , there is an 
ambiguity in carrying out the series expansion. We could either use the form given above 
and expand the denominator in a series expansion in e~ 2 ™ so that the criterion described 
above is satisfied, or express it as e 2mv / (1 — e 2mv ) 2 and expand it in a series expansion in 
e 2mv in which case the sum over p will be bounded from below rather than from above. 
One finds that depending on the angle between S 1 and S 1 , only one of these expansions 
produces the degeneracy correctly via (1 2 . 2 2 [) [1]. The physical spectrum actually changes 
as this angle passes through 90° since at this point the system is only marginally stable. 
On the other hand our degeneracy formula (12.141) . (12. 15ft implicitly assumes that we have 
expanded this factor in powers of e~ 2mv since only in this case the sum over p in ( 12.221) is 
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bounded from above for fixed m, n. Thus as it stands the expression for d(Q, P) given in 
(12.141) . (12.151) is valid for a specific range of values of the angle between S 1 and S , which, 
in the second description of the system, corresponds to the sign of the axion field. For 
the other sign of the axion we need to carry out the integral over a different contour C 
defined as 

p 2 = M u a 2 = M 2 , v 2 = M 3 , 
< pi < 1, < ai < N, < vi < 1 , (2.29) 

where Mi, M 2 and M 3 are large positive numbers with M 3 << Mi, M 2 . In the convention 
that we shall be using in this paper, we need to use the contour C for positive sign of the 
axion and the contour C for negative sign of the axion. 

It turns out that walls of marginal stability, - codimension one subspaces of the asymp- 
totic moduli space on which the BPS mass of the system becomes equal to the sum of 
masses of two or more other states carrying the same total charge, - are quite generic 
for quarter BPS states in M = 4 supersymmetric string theories [11], and we expect the 
spectrum to change discontinuously as the asymptotic moduli fields pass through any 
of the walls of marginal stability^ Thus the expression for the degeneracy given in this 
section holds only in a finite region of the moduli space, bounded by the walls of marginal 
stability. This will be discussed in more detail in £j3l 

Another point about the formula (I2.14p is that although it was derived for special 
charge vectors Q, P described in (12.121) . it has been expressed as a function of the T- 
duality invariant combinations P 2 , Q 2 and Q ■ P. Even though we expect T-duality to be 
a symmetry of the theory, it is not guaranteed that the formula written in this fashion 
hold for all charge vectors. First of all a T-duality transformation acts not only on the 
charges but also on the asymptotic moduli. Had the spectrum been independent of the 
asymptotic moduli, we could have demanded that the spectrum remains invariant under 
T-duality transformation of the charges. However if a T-duality transformation takes 
the asymptotic moduli fields across a wall of marginal stability, then all we can say is 
that the spectrum remains unchanged under a simultaneous T-duality transformation of 
the moduli fields and the charges, but if we are sitting at a fixed point in the moduli 
space then the spectrum is not invariant under T-duality transformation on the charges. 

4 If a state becomes marginally stable on a surface of codimension > 2, then we can always move around 
this subspace in going from one point to another and hence the spectrum cannot change discontinuously. 
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Second point is that even if we ignore the issues related to the walls of marginal stability, 
two charge vectors carrying the same values of P 2 , Q 2 and Q ■ P may not necessarily be 
related by a T-duality transformation^ In that case the degeneracy of states for these two 
charge vectors could be different. An example of this is that a state that carries fractional 
momentum along S l can never be related to a state carrying integer momentum along 
S 1 , although they may carry same values of Q 2 , P 2 and Q ■ P. Both these issues, together 
with the S-duality transformation properties of the degeneracy formula, will be discussed 
in gl 

3 Walls of Marginal Stability 

As has been briefly mentioned in the degeneracy formula given in (12.141) . (12.151) is 
expected to be valid within a certain region of the moduli space bounded by codimension 
one subspaces on which the BPS state under consideration becomes marginally stable. As 
we cross this subspace of the moduli space, the spectrum can change discontinuously. In 
this section we shall study in some detail the locations of these walls of marginal stability 
so that we can identify the region within which our degeneracy formula will remain valid. 

Let us consider a state carrying electric charge Q and magnetic charge P and examine 
under what condition it can decay into a pair of half-BPS states. This happens when its 
mass is equal to the sum of the masses of a pair of half BPS states whose electric and 
magnetic charges add up to Q and P respectively. Since for half BPS states the electric 
and magnetic charges must be parallel, these pair of states must have charge vectors of 
the form (aM, cM) and (bN, dN) for some constants a, b, c, d and a pair of r-dimensional 
vectors M, N. We shall normalize M, N such that 

ad — be = 1 . (3-1) 

Then the requirement that the charges add up to (Q, P) gives 

M = dQ-bP, N = -cQ + aP. (3.2) 

Thus the charges of the decay products are given by 

{adQ - abP,cdQ - cbP) and (-bcQ + abP, -cdQ + adP) . (3.3) 

5 Generically they are related by a continuous T-duality transformation but only a discrete subgroup 
of this is a genuine symmetry of the theory. 
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Note that under the scale transformation 

a b\ fa b\ ( X 



c d \c d \0 X- 1 



(3.4) 



eqs. (l3.ip and (13.31) remain unchanged. There is another discrete transformation 

(3.5) 



a b\ fa b W 1 
c d) ~* \c d) \-l 

which leaves (13 .ip unchanged and exchanges the two decay products in (13.31) . A pair of 

matrices ( ^ ^ j related by (13.41) or (13. 5p describe identical decay channels. 

In order that the charge vectors of the decay products given in (13.31) satisfy the charge 
quantization rules we must ensure that aM = adQ — abP and bN = —bcQ + abP belong 
to the lattice of electric charges and that cM = cdQ — cbP and dN = —cdQ + adP belong 
to the lattice of magnetic charges. For the charge vectors Q, P given in (12.121) this would 
require 

ad, ab, be e 71, cd e N 7L . (3.6) 

The condition cd G N 7L comes from the requirement that cdQ — cbP is an allowed 
magnetic charge. In particular for a Q of the form (12.121) . a magnetic charge cdQ represents 
a state with Kaluza-Klein monopole charge —cd associated with S 1 . Since this charge is 
quantized in units of N, cd must be a multiple of N. We shall denote by A the set of 

matrices f ° subject to the equivalence relations (13.41) . (13.51) and satisfying (13. ip . 

ASH). 

It is instructive to determine the structure of the set A. We shall first show that 
using the scale transformation (13. 4p we can always choose a, b, c and d to be integers and 
furthermore the solution is unique for given ad, ab, be and cd. Since ad, ab, be and cd are 
all integers, we can express them as products of prime factors: 

Wb\ = l[p?, \cd\ = l[pt', \ad\=Y[ P r, \bc\ = Y[ P r, (3.7) 

i i i i 

where the product over i runs over the prime numbers pi and Ti, Sj, Ui and Vi are non- 
negative integers satisfying 

r i + s i = u i + v i Mi. (3.8) 

Furthermore since ad and be differ by 1, they cannot have a common factor. This shows 
that either Ui or V{ must vanish: 

U{ v i = Wi . (3.9) 
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Let us now look for integer a, b, c and d satisfying (13. 7p . For this we use the ansatz: 

\a\ = U P T, N = IK 1 ' l c l = n^' \d\ = Hpt, (3-10) 

i i i i 

where a i: bi, Ci and di are non-negative integers. Eq. (13.71) now gives 

ai + bi=ri, d + di = Si, a { + d { = m, b i J rc i = v i . (3.11) 

Now f)3.9p tells us that for any given i either Ui or Vj (or both) are zero. If Ui — then 
the only possible solution to (13. lip is 

a, = 0, di = 0, bi = r h Q = . (3.12) 

On the other hand if v , = then we must have 

bi = 0, Ci = 0, cii = n, di = Si. (3.13) 

This gives a unique expression for \a\, \b\, \c\ and \d\ using prime factorization. Up to an 
overall factor of —1 which can be removed with the help of the scale transformation (13. 4p . 
we can determine the signs of a, b, c and d in terms of the signs of ab, cd, ad and be. 

The requirement that cd is a multiple of N implies that the prime factors of iV are 
shared by c and d. In case N is prime either c or d must be a multiple of N. Using 
the freedom (I3.5P we can ensure that c is a multiple of N. In this case the matrices 

f c describe elements of r (iV) modulo multiplication by — lj*j Since T (2) = 1^(2), 

for N = 2 the set of matrices ( a ^ j may be identified as the elements of Fi (N) modulo 
multiplication by —1. On the other hand using the freedom of multiplication by —1 we 
can convert any T (3) matrix to a 1^(3) matrix. Thus for N = 3 the matrices 

may be chosen to be Ti(3) matrices. 

The case N — 1 with M. = K3, corresponding to heterotic string theory on T 6 in the 
second description, is somewhat special. In this case the set A consists of PSL(2, 2L) 
matrices subject to the equivalence relation (13. 5p . 

We shall now determine the wall of marginal stability corresponding to the decay 
channel given in (13. 3p . Our starting point will be the formula for the mass m(Q,P) of a 



6 To(A^) contains matrices ( ° , ) satisfying ad — be = 1, a,b,d S Z, c £ JV 
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BPS state carrying electric charge Q and magnetic charge P [20,21] 

1 



m(Q,Pf = -^-(Q-fooPfiM^ + L^Q-r^P) 

>->oc 



+2 [(Q T (M oc + L)Q)(P T (M oc + L)P) - (P T (M oc + L)Q) 2 ] 1/2 , 

(3.14) 

where r = a + iS and the subscript oo denotes asymptotic values of various fields. This 
expression is manifestly invariant under the T- and S-duality transformations described in 
eqs. fl2.5p - fl2.7p and fl2.9p - fl2.lip . In order that the state (Q, P) is marginally stable against 
decay into (adQ — abP, cdQ — cbP) and (—bcQ + abP, —cdQ + adP), we need 

m(Q, P) = m{adQ - abP, cdQ - cbP) + m(-bcQ + abP, -cdQ + adP) . (3.15) 

Using ( 13.141) . ( 13.15P and some tedious algebra, we arrive at the condition 

(■ °l 6c ) 2 +k + Ay^ (1+ ^, (3.i6) 

where 



2cd J V 2cdJ Ac 2 d? 



E = cd(Q T (M oc + L)Q) + ab(P T (M OQ + L)P) - (ad + 6c)(P r (M oc + L)Q) 
[(Q T (M 00 + L)g)(P^(M 00 + L)P) - (P T (M oc + L)Qf} 1/2 

Note that E depends on M^, the constants a,b,c,d and the charges Q, P, but is inde- 
pendent of T^. Thus for fixed P, Q and M^, the wall of marginal stability describes a 
circle in the (a^, S^) plane with radius 



R= vT+E*/2\cd\, (3.18 

and center at 

C= (- 

2cd ' 2cd 
This circle intersects the real at 

a/c and b/d. (3.20) 

The cases where either c or d vanish require special attention. First consider the case 
c = 0. In this case the condition ad — bc = l implies that a — d — 1. By taking the c — > 
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limit of (13.161) . (I3.17P we see that the wall of marginal stability becomes a straight line in 
the (aoo, Soo) plane for a fixed M^: 



b(P T (M QO + L)P) - (P T (M 00 + L)Q) 



[(g T (M 00 + L)g)(P r (M 00 + L)P) - (P T (M oc + L)Q) 



21 1/2 



S n 



0. (3.21] 



On the other hand for d = we have be = — 1 and we can choose b — —1, c — 1. In the 
cf — > limit of (I3.16p . (13.171) we get another straight line 

a(P T (M 00 + L)P)-(P T (M 



L)Q) 

[(Q T (M 00 + L^P^M^ + L)P) - {P T {M 00 + L)Q) 



21 1/2 



5*^ 



0. (3.22) 



We now notice that this has exactly the same form as (13.211) with b replaced by a. Thus 

these do not give rise to new walls of marginal stability. In fact the c = and d = cases 

are related by the equivalence relation (13. 5p . 

In order to get some insight into the geometric structure of the domain bounded by 

these marginal stability walls it will be useful to study the possible intersection points of 

these walls in the upper half plane. Let us consider a pair of such walls characterized 
ai bi 



by the matrices 



Ci di 



and 



a 2 b 2 
c 2 d 2 



. A convenient procedure for studying their 



intersection is to convert one of them (say the first one) to a straight line by an SL(2, 
transformation. We define 



diT — bi 



(3.23) 



— C\T + ai 

Then it is easy to see that in the (a'^, S^J plane the two walls get mapped to the curves 

P T (M OQ + L)Q 



a' + 



(Q T (M oc + L)Q)(P r (M 00 + L)P) - (P r (M 00 + L)Qf 



Q> 

1/2 oo 







(3.24) 



and 



where 



ad + be 
2cd 



+ SL 




A^d 2 



1+P 2 



Q = d x Q-b x P, P = -c 1 Q + a 1 P, 



d\ 



-bi 
a i 



a 2 b 2 
c 2 d 2 



(3.25) 

(3.26) 
(3.27) 
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; , _ cd(Q T (M oa + L)Q) + ab(P T (M oc + L)P) - (ad + bZ)(P T (M QO + L)Q) 

h = TZ Z~ZZ Z Z 11/2 ' W-^) 

(Q T (M OQ + L)Q)(P^(M OQ + L)P) - (P T (M 00 + L)Qf 

If either c or d vanishes i.e. if either a\jc\ = a 2 /c 2 or a\jc\ = b 2 /d 2 then (I3.25P reduces 
to a straight line of the form (13.211) or (13 .22 j) with P, Q, b (or a) replaced by P, Q, b (or 
a) respectively, and one finds that the only point of intersection of this line with (I3.24p in 
the upper half plane is at zoo. In the r^, plane this corresponds to the point a±/ci. If 
neither c nor d vanishes, then by eliminating a'^ from (13.241) and ( 13.251) we get 

, c , , 2 Q T (M 00 + L)QP T (M 00 + L)P 

Wo 



'00/ 



(Q T (M 00 + L)Q)(P T (M OQ + L)P) - (P T (M QO + L)QY 

g^oo + L)Q + (^/cd)P T (M 00 + L)P 

(Q^M^ + L)Q)(P T (M 00 + L)P) - (P T (M 0O + L)Q)A ^ 



(5a + ^-i =0< 

&d 2 

Using the conditions ad — bc = 1 and a,b,c,d€L 7L we see that ad and be have same signs 
and also that (ad + be) 2 > 1. As a result for 5^ > 0, each term in the left hand side of 
(13.291) is non-negative and the only possible solution is 

6 = or a = , 5^ = 0. (3.30) 

Going back to the original variables using (13.231) . (13.271) we see that these correspond to 
the following two cases: 

h «2 h b 2 h 

-r = — or ~r = -r, = — . 3.31 

d\ c 2 d\ d 2 d\ 

Collecting all the cases together we see that the pair of circles characterized by the matrices 
' 1 and \ 2 ^ 2 ) never intersect in the interior of the upper half plane and 



ci d\ J \c 2 d 2 

intersect on the real axis if and only if the sets 



and (3.32) 



ci d\ J [ c 2 d 2 

have an overlap. 

From this analysis we see that the set of marginal stability walls divides the upper half 
Too plane into many domains, with each domain bounded by a set of walls described by a 
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set o f _ <C: b J)>- S.,oe the wa„ S m eet on the tea, axis ot at »co, the regions 



have possible vertices on the real axis or at zoo, but never in the interior of the upper half 
Too plane. We also note that while the shape of the walls in the Too plane depends on the 
charges and the moduli M m , their intersetion points are determined only in terms of the 

associated matrices ( a% ^ 

\Ci di 

How do we determine the walls bordering a given domain? We shall illustrate this 
for the domain that includes the large region and is bounded on the left (right) by 
the line ( 13.21 j) for b = (6 = 1). Since these lines intersect the real axis at and 
1 respectively, the domain is bounded from below by set of circle segments in the upper 
half plane described by (13.161) . of which the first one begins at and ends at some point 
Pi on the positive real axis, the second one begins at Pi and ends at some point P2 to 
the right of Pi etc. with the final segment ending at 1. The parameters a, b, c, d for each 
segment must be chosen such that given the starting point, it travels maximum possible 
distance to the right; any other segment that travels less distance will lie beneath the 
maximal segment and will not be of relevance for this computation. Thus in order to 
determine these segments we first need to study the following general question: given a 
rational number p/q > 0, what is the maximum distance a wall can travel given that it 
starts at p/q! 

Let us first consider the case p/q > 0. We choose p, q to be relatively prime and p > 0, 
q > 0. Now eq. (l3.20p shows that the circle (13.161) intersects the real axis at a/c and b/d. 
Thus either a/c or b/d should be equal to p/q. Without any loss of generality we can 
choose a/c = p/q and hence a = p, c = q up to an irrelevant overall sign. The constraint 
ad — 6c = 1 now gives 

b p 1 , 

1 = A ■ 3 - 33 

d q qd 

We need to choose d so as to maximize the right hand side. Thus d must be chosen to 
be the maximum possible negative number subject to the requirement that cd = qd is a 
multiple of N and b = (pd — l)/g is an integer. 

In the special case where p/q = 0, we have a/c = and hence a = 0. The condition 
ad — be = 1 tells us that we have (up to an overall sign) b = — 1, c = 1. The second 
intersection point b/d = — 1/d is maximized for maximum negative value of d subject to 
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1 1/2 1 1/3 1/2 2/3 1 

N=l N=2 N=3 



Figure 1: A schematic diagram representing the domain 1Z in the upper half r M plane, 
bounded by the walls of marginal stability, for ZSj., ^ an d ^3 orbifolds. The shapes 
of the circles and the slopes of the straight lines bordering the domain depend on the 
charges and other asymptotic moduli, but the vertices are universal. 

the requirement that cd is a multiple of N. This gives d = —N and hence 

Following this rule we can now construct the circles which border the domain from 
below. For the purpose of illustration we shall carry out the first few steps. Using (13.341) 
we see that the first circle segment beginning at ends at 1/N. For N = 1 this completes 
the story since we have already reached the point 1. For N > 2 we need to proceed further. 
Taking p/q — 1/N we see from ( 13 .33|) that b = (d— l)/N. Since q = N, the condition that 
qd is a multiple of N is trivially satisfied. Thus we need to choose d to be the maximum 
negative number for which b = (d — 1)/N is integer. This gives d = —(N — 1) and b = — 1. 
Thus the circle ends at b/d = 1/(N — 1). For N = 2 this completes the story but for 
iV > 3 we need to proceed further. At the next stage we begin with p/q = 1/ (N — 1) and 
get b = —(N - 1), d = -N(N - 2). Thus b/d = (N — 1)/((N(N - 2)). This does not 
complete the story for any N > 3; e.g. for iV = 3 this gives b/d = 2/3. By continuing this 
process one can show that for N = 3 we reach the point 1 at the next step via the wall 
corresponding to ( ° \ J = ( \ \ \ but for higher iV the story continues further. A 



c d I V 3 —1 

schematic diagram representing these domains for N = 1, 2 and 3 have been shown in 
Fig. [H As will be discussed in £|6l for N > 4 the number of such walls is infinite. 

This finishes our general analysis of marginal stability walls and domains bounded 
by them. Now we focus on one particular domain, - the one in which the degeneracy 
formula given in (I2.14p . (12. 15[) is valid. Since the calculation leading to (12.141) . (12.151) 
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was performed in a specific corner of the moduli space, - the weakly coupled type IIB 
string theory, - all we need to know is how this region is situated with respect to the 
various marginal stability walls described here. However to address this issue we need 
to determine the relation beween the moduli parameters in the first description, - as an 
orbifold of type IIB string theory, - and the moduli t^, M^, appearing in the BPS mass 
formula (13.141) . We have already stated that r M = a 0O + iS' 0O denotes the asymptotic value 
of the complex structure of the torus (S 1 x S 1 )/ TL^. We have worked in a region of the 
moduli space where it is finite, i.e. is neither too large nor too small. The relation 
between the other moduli fields in the first description and the matrix valued moduli 
field M can be found by following carefully the duality chain that takes the theory to its 
second description, and using the identification of M with the geometric quantities in this 
description as given in (12.31) . Let us denote, in the first description of the theory, by g 
the ten dimensional coupling constant, by {2 , K\fa') A V the volume of A4 measured in the 
string metric, and by (2n\/a') 2 A the area of S 1 x S 1 / %n measured in the string metric. 
Let us also set the anti-symmetric tensor field and all the Ramond-Ramond fields to zero. 
Then using the standard duality transformation rules one can show that the relevant 4x4 
component part of the matrix M that couples to Q and P are given by 

(l/V \ 



Moo 



(3.35) 



g*/(A*V) 

V 

V A 2 V/g 2 J 
Using the expression for Q and P given in (12.121) one now gets 

w - &®f + f 

P^M^ + ^P = iQl ~° 2 + ^ V J 2 + V + 2(Q 1 -C), 

P T (M 00 + L)Q = -£~ + J- (3-36) 

The computation of the degeneracy was done in a region where g is small but A and V 
are of order 1. Using (13.361) one can now see that in this region 

Q T (M OQ + L)Q » P T (M OQ + L)P, |P T (M 00 + L)Q\ . (3.37) 

Even when we deform M m away a little from its diagonal form (I3.35f) . eq.f l3.37n continues 
to hold. 



23 



We now study the implication of (13.371) on the location of the walls of marginal stability 
in the Too plane. First consider the case where d ^ 0. In this case we see from ( 13. 17ft . 
( 13371) that 

~d-y PT( Moc + L)P >>X - (3 ' 38) 
Hence for the circle in the plane described by (13.161) we have 



^ <} iSr-£i <<l - (3 - 39) 

In other words in the (doo, Soo) plane the uppermost point on the circle (13.161) lies little 
above the = axis and its center lies deep down in the lower half plane. Since ~ 1 
in the region of the moduli space in which we have worked, we see that this region lies 
above all the circles described by ( 13. 16[) in the plane. 

Next we consider the case c = 0, - as discussed earlier the d = case is equivalent to 
this. In this case eq. ( 13. 36ft shows that the coefficient of Soo in ( 13.211) is small for finite b 
and small g. Thus for fixed M^, these lines are almost vertical in the (a^, Soo) plane and 
are given by 

(3.40) 

The 6 = line corresponds to the wall of marginal stability described at the end of §2J 
and as we described in £j2j the degeneracy formula actually jumps across this wallQ Since 
our degeneracy formula has been derived in the small region, we see that the region of 
validity of our formula is bounded by the b = — 1 line on the left, 6 = 1 line on the right 
and a set of circle segments below. 

To summarize, the region of the moduli space in which we have carried out our analysis 

consists of two domains. One of them, lying between the ( \ J = ( \ \ \ line and 



c rf) = (o 1 ) ^ ne ex ^ enc ^ s ^° large region, and is bounded from below by 

a set of circles. For later reference we shall call this the right domain 1Z and denote by Br 

a b 



the set of matrices ( ^ ] labelling the boundaries of this domain. The other domain, 



7 In the analysis of [1] the jump occured exactly across the line doo = 0. This can be traced to the fact 
that there the analysis was carried out using weak coupling approximation and did not take into account 
the backreaction due to switching on J. Indeed for J = we have P T {Moo + L)Q = 0, and (13.21|) for 
b = reduces to the vertical line aoc = 0. 
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between th e( : *) - g ?) d £ »)- (J J ) line a.so extends 

to the large region and is bounded from below by a set of circles. We shall call this 
the left domain £, and denote by Bl the set of matrices ( a \ ] labelling the boundaries 



c 

of this domain. We shall argue in §[5] that the set Bl is obtained simply by multiplying 

1_1 



the elements of Br by the matrix I ^ ^ J from the left. 

One important issue that we would like to address is: how does the degeneracy formula 
change as we move across a wall of marginal stability? We already know that as we move 
across the wall (I3.2ip for b = the degeneracy formula continues to be given by (12.141) 
except for a change in the location of the integration contour from C to C. We shall argue 
in §H] that this is a general phenomenon; as we move across any line of marginal stability 
the degeneracy formula will be given by the same expression except for a change in the 
integration contour. 

4 Duality Transformation of the Degeneracy Formula 

As noted in §0 the degeneracy formula (12.141) . (I2.15P has been written in terms of T- 
duality invariant combinations Q 2 , P 2 and Q ■ P although we have derived the formula 
only for a special class of charge vectors. In this section we shall discuss what information 
about the degeneracy formula can be extracted using the T- and S-duality symmetries of 
the theory. 

We begin by studying the consequences of the T-duality symmetries of the theory. It 
follows from (12.51) . (I2.6P that if a T-duality transformation takes a charge vector (Q,P) 
to (Q', P') then 

Q' 2 = Q 2 , P' 2 = P 2 , Q'-P' = Q-P. (4.1) 

However there may be pairs of charge vectors with the same Q 2 , P 2 and Q ■ P which 
are not related by a T-duality transformation. Clearly T-duality invariance of the theory 
cannot give us any relation between the degeneracies associated with such a pair of charge 
vectors. In what follows we shall focus on charge vectors (Q', P') which are in the same 
T-duality orbit of a charge vector (Q, P) for which we have derived (I2.14p . 

We have denoted by TZ the right domain of the region of the moduli space described 
in £j3]in which the original formula for d(Q, P) is valid. It is bounded by a set of marginal 
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(X b 

stability walls labelled by ( , ) G Br. Let TZ' denote the image of TZ under the T 



c d 

duality map. In this case we expect d(Q', P') in the region TZ' to be equal to d(Q, P) 
given in (12.141) : 

d(Qi,P>) = - [ dpdGdve-^ N ~ pQ2+ ° p2 l N+2 ^ „ 4r- > 

= I / dp do dve~™ (NpQ'^P^/N + 2vQ'. n (4.2) 

ic $(p, a, 

where C has been defined in ( 12.151) . In going from the first to the second line of (14.21) we 
have used (14. ip . 

Let us now determine the region TZ' . Since under a T-duality transformation M — > 
QMQ T , and since TZ' is the image of 71 under this map, 71' is bounded by walls of marginal 
stability described in (13.161) . (13.171) with f ° G Br and in (13.171) replaced by 

n~ l M 00 (pt r )- 1 . Using §ZM) we see that this effectively replaces (Q, P) by (Q', P') in 
(13.171) . Thus 7Z' is the region of the upper half plane bounded by the circles: 



2cd J \°° 2cdJ 4c 2 d 2K ; ' \c d 

where 

^ cdjQ'TjM^ + L)Q') + ab{P^iM^ + L)P') - {ad + &c)(P^(M QO + L)QQ 
[(Q ,T (M 0O + L)g')(^ /T (M 00 + L)P') - (P' r (M 00 + L)Q') 2 ] 1/2 

(I4.2p - (I4.4I) are valid for any charge vector (Q',P') which can be related to the charge 
vectors given in (I2.12p via a T-duality transformation. 

Next we shall analyze the consequences of S-duality symmetry. An S-duality transfor- 
mation changes the vector (Q, P) to another vector (Q", P") and r to r" via the formulae 
( 12. 9p . (12.101) . Thus if TZ" denotes the image of the region TZ under the map (I2.10p . then 
S-duality invariance implies that inside TZ" the degeneracy d(Q", P") is given by the same 
expression (12.141) for d(Q,P): 

d ($» } p") = 1 f d p d?rd y e -™(NpQ^PyN + 2vQ-P) __L_ . (4.5) 
N Jc &(p,cr,v) 

We would like to express the right hand side of (14.51) in terms of the vectors Q" and P" . 
For this we define 

(; s)-U t)*™- (4 - 6) 
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and 

a 2 (3 2 -2a(3 



p 










l-l 


a'[ + io'{ 


H 












7 2 5 2 __ z 2^5 | | a | . (4.7) 
—67 — /55 (a5 + P'y) 

Using (12. 9p . (14. 6p . (14.71) one can easily verify that 

e -ni{NpQ 2 +aP 2 /N+2vQ-P) _ e -7ri(7Vp"Q" 2 +3 : "P" 2 /A r +2y Q"-P") U 

and 

dp da dv = dp" da" dv" . (4-9) 
Furthermore, with the help of eq. (12.271) one can show that [3] 

$(f?',a",v") = $(p,a,v). (4.10) 

If C" denotes the image of C under the map (14.71) then eqs. (14.8(1 - (14. lOj) allow us to express 
(143(1 as 

d (Q» p") = 1 / d 7/> fa" fa" e -MN?>Q''*+Z''P''VN + ZvQ».P») ^ 1 (4 n) 

N Jc» $(p",a",v") 

To find the location of C" we note that under the map (14.71) the real parts of p, a and v 
mix among themselves and the imaginary parts of p, a and v mix among themselves. The 
initial contour C corresponded to a unit cell of the cubic lattice in the (pi,cri,vi) space 
spanned by the basis vectors (1,0,0), (0, N, 0) and (0,0,1). The unimodular map (14.71) 
transforms this into a different unit cell of the same lattice. We can now use the shift 
symmetries 

$(p,a,v) = $(p + l,a,v) = $(p,a + N,v) = $(p,a,v + l) , (4.12) 

which are manifest from (12. 16|) . to bring the integration region back to the original unit 
cell. Thus C" and C differ only in the values of the imaginary parts of p, a and v. Using 
(12.151) . (14.711 we see that for the contour C", 

Pz = a 2 M 1 + p 2 M 2 + 2a(3 M 3 , 
a? % =fM 1 + 5 2 M 2 + 275 M 3 , 
v » = -ary M x -(3~8M 2 - (ad + ^7) M 3 . (4.13) 

Thus C" is not identical to C, - a fact first noticed in [13]. We could try to deform C" 
back to C, but in that process we might pick up contribution from the residues at the 
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poles of <&(p" ,a",v"). Thus we see that the degeneracy formula (14. lip for d(Q",P") is 
not obtained by simply replacing (Q,P) by (Q",P") in the expression for d(Q,P). The 
integration contour C also gets deformed to a new contour C" . 

Let us now analyze the region TZ" of the asymptotic moduli space in which (14. lip is 
valid. This is obtained by taking the image of the region 71 under the transformation 
(12.101) . To determine this region we need to first study the images of the curves described 
in (13.161) in the — Soo plane. A straightforward analysis shows that the image of (13.161) 
is described by the curve 



a"d" + b"c"\ 2 / E" V 1 M 



where 



(4.15) 



a" b"\ _ (a (3\ (a b 
c" d") ~ V7 SJ \c d 

and 

„ _ c"rf // (g /,T (M 00 + L)Q") + a"V'{P" T (M 00 + L)P") - {a"d" + b"c")(P" T (M QO + L)Q") 



E 



[(Q"' r (M 00 + L)Q")0P" T (il4o + L)P") - (P" r (M 00 + L)Q") 2 ] 1/2 

(4.16) 

( a" b" \ 

It can be easily seen that I ^ ^, J satisfy the relations (13. ip . (13.61) . and the equivalence 

/ a " y \ 

relations ( 13 .4h . f]3. 51) translate to identical equivalence relations on I ^ ^„ J. Thus the 

/ a " y 

collection of all the matrices I „ ,„ ) describes same set A as the collection of all the 

matrices ( ° \ 
c d 

Now recall that the original region TZ was bounded by a set of marginal stability walls 
c J ^ r ^^ ms ^ e re §i° n ^ s bounded by the collection of walls described 

by (33U), with d") G (" 

At this stage it will be instructive to compare the expressions for the degeneracies 
d(Q,P) and d(Q",P"). There are two key differences. First of all although both are 
represented by similar looking integrals, the contours of integration are different in the 
two cases. Second the region of validity TZ" of the expression for d(Q, P) is not obtained by 
simple replacement of (Q, P) — > (Q", P") in the expression for d(Q, P); the original region 
of validity of our formula was bounded by a set of marginal stability walls corresponding 
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to a set of matrices G Bit, whereas the new region of validity of the formula is 

bounded by a set of walls corresponding to the matrices \ „ \„ ) G ( a ^ J £>_r. 



Typically the new charge vectors (Q", P") are not of the type (I2.12p . but in some cases 
it is possible to obtain the charge vectors (Q", P") by T-duality transformation of another 
charge vector (Q, P) of the form given in (12.121) satisfying 

Q 2 = Q" 2 , ?* = P" !i , Q.p = Q ».p» (417) 

Then for d(Q,P) our original formula for degeneracy holds. Using the results given in 
L2]) - (14.41) we can now conclude that in the region TZ" bounded by the walls (I4.16P with 



a" b" 
c" d" 



eB R , (4.18) 



the degeneracy d(Q",P") is given by (14. lip with the contour C" replaced by C. 

Thus we now have expressions for d(Q", P") in two different domains, - TZ" and TZ". 
In both regions the degeneracy is given by an integral. The integrand in both cases are 
same, but in one case the integration contour is C while in the other case it is C" . This 
shows that as we cross the walls of marginal stability to move from the region TZ" to TZ" 
the expression for d(Q", P") changes by a modification in the location of the contour of 
integration. 

Even though the result was derived under the assumption that (Q", P") can be related 
to a charge vector of the type fl 2 . 1 2 [) by T-duality, it is natural to assume that this 
phenomenon is more general. In particular, a natural postulate will be that as we move 
from the domain TZ corresponding to the set Br to another domain corresponding to the 
set Br, with G Tx(N), the expression for d(Q , P) gets modified by a 

replacement of the integration contour C to the contour C" corresponding to 

p 2 = a 2 m 1 + p 2 m 2 + 2ap m 3 , 

a 2 = 7 2 Mi + S 2 M 2 + 2j5M 3 , 
V2 = -5q Mi - j35M 2 - (&S + fly) M 3 
< pi < 1, < G X < N, < vx < 1 , (4.19) 

Mi, M 2 and M3 being large but fixed positive numbers with M3 << Mi,M 2 and a = 5, 
]3 = y/N, y = N(3,5 = a. 
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Can every domain be related to TZ or C via S-duality transformation? If this is so then 
we can use duality invariance together with the information about the contours C and C 
appropriate for TZ and C, - as given in (12.151) and (12.291) . - to find out the location of the 
integration contour in every other domain. For N = 1, 2, 3 we can answer this question 

in the affirmative as follows. In these cases the matrix ^ a ^ J labelling any wall can 
be represented by an element of Ti(N) and hence can be related to the matrix 
via a duality transformation of the form ( 14.151) . Since describes the boundary 

between the domains TZ and £, the two domains bordering the wall ( , J must be 



c d 

related to the domains TZ and £ by the same duality transformation. Since this can be 
done for every wall of marginal stability, we see that every domain must be related to 
either TZ or C via a duality transformation. 

For N > 4 not all the walls can be related to the wall ( ^ ^ ) via duality transfor- 



x x , 

mation. Nevertheless it may still be possible to relate all the domains to TZ and L via 
duality transfomation, with the walls not related to ( _ \ \ being related to some other 



x x , 

wall of TZ or C We shall not attempt to settle this issue here. 

Given two contours on two sides of a marginal stability wall, if we try to deform one 
to the other then typically we shall encounter poles of the integrand, and the two results 
will differ by the sum of the residues at these poles. To get a physical insight into how 
much this change is as we move across a given wall, let us consider the wall corresponding 

to the matrix ( !? J . This is the wall separating the regions TZ and C, and across this 



wall the original dyon becomes unstable against decay into states with charges (Q, 0) and 
(0, P). Since across this wall the contour changes from constant negative t>2 to constant 
positive t>2, we pick up the residue at the pole v = 0. Now near v = the function $ has 
the behaviour [2,3,5,6] 

&(fi,a,v) ^ —A^v 2 f(Np)g(a/N) + 0(v 4 ) , (4.20) 

where f(r) and g(r) are two functions which have the interpretation of inverse of partition 
functions associated with electric and magnetic half BPS states of the theory. Performing 
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the integral over v in f)2.14p around the v — point now gives 

- (Q ■ P) { jf * dpe-™ N ^ 2 (/(iVp))- 1 } |1 £ dae-^ p2 l N (g(a/N))-^ . (4.21) 

This formula can be given a simple physical interpretation. The second and the third 
factors represent the degeneracies of the electric and magnetic half BPS states into which 
the original dyon decays on the marginal stability wall. The Q ■ P factor on the other 
hand is associated with the supersymmetric quantum mechanics describing the relative 
motion of the electric and the magnetic system [1] . In particular it represents the number 
of states whose binding energy vanishes as we reach the marginal stability wall [22,23]. 
Thus this is the number of states which disappear from the spectrum as we cross the wall. 

5 The Large Charge Limit 

In this section we shall argue that even though the complete spectrum changes discon- 
tinuously when the asymptotic value of the axion field changes sign, the large charge 
expansion is not affected by this change. Our starting point is the result derived in [1-3] 
that the poles in the integrand in the expression for d(Q, P) come from the second order 
zeroes of $ at§ 

n 2 (ap — v 2 ) + jv + n{a — m{p + m 2 = , 

f 1 
mini + m 2 n 2 + — = - , 

mi e NIL, ni,m 2 ,n 2 e^, j'621 + 1. (5.1) 

For large charges the leading contribution to the degeneracy comes from the pole at 
n 2 = 1 [1-5,24]. Contribution from the poles with n 2 ^ 1 are exponentially suppressed 
compared to the leading contribution. 

Consider now the contours C and C" given in (I2.15P and (I4.19p respectively. Both 
contours have the same range of integration over the real parts of p, a and v 

< pi < 1, < a x < N, 0<v<l. (5.2) 

8 For M. = K3 this result can be found in appendix E of [1] . For the general case the set of all the 
zeroes and poles of $ were listed in [3], but we did not attempt to separate out the zeroes from the poles. 
A careful analysis shows that the only zeroes come from the set (|5.ip . 
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However they differ in the values of p 2 , u 2 and v 2 . We have for 



C : p 2 = Mi, a 2 = M 2 , v 2 = -M 3 , 

C" : p 2 = a 2 M l + (3 2 M 2 + 2a(3 M 3} a 2 = f M 1 + 5 2 M 2 + 2^5 M 
v 2 = -dry M l - f35M 2 - (a5 + ^7) M 3 , 



(5.3) 



where Mi, M 2 and M 3 are fixed numbers with M 1} M 2 » 1 and |M 3 | << M 1; M 2 . Our 
goal will be to show that we can deform the contour C to C" without hitting any of the 
poles given in (15.11) except those with n 2 = 0. Since for large charges the contribution 
to the degeneracy comes from the poles at n 2 = 1 up to exponentially suppressed terms, 
this will show that the change in the degeneracy due to the change in the contour of 
integration is exponentially suppressed compared to the leading contribution. 

To proceed we shall choose a specific path along which we deform the contour. We take 
this to be along the straight line joining the points (M l5 M 2 , M 3 ) and (a 2 Mi + (3 2 M 2 + 
2af3 M 3 , 7 2 M x + 5 2 M 2 + 2^5 M 3 , -Sty M x -p8M 2 - (a5 + f3^) M 3 ) in the (p 2 , a 2 , v 2 ) space. 
Points on this line can be parametrized by a real number A lying between and 1, with 



For each A we take the real parts of p, a and v to be in the range (15.21) . 

We want to show that (15.11) and (15.41) have no simultaneous solutions except for n 2 = 0. 
For this we write down separately the real and imaginary parts of eq. (15.1 j) . as well as the 
constraint on m, n, j: 



Mi + A {(a 2 - 1) M x + (3 2 M 2 + 2S(3 M 3 } 

M 2 + A {7 2 Mi + (5 2 - 1) M 2 + 275 M 3 } 

-M 3 + A {-«7 Mi -pSM 2 - (&8 + ^7 - 1) M 3 } . 



(5.4) 




(5.5) 




mini + m 2 n 2 + — = - . 



(5.6) 



(5.7) 



We now eliminate mi and m 2 using (15.61) . (15. 7h to write (15.51) as 



(a 2 p 2 -v 2 )+ (j ni) A 



( 



J 

m 



) 



+ aj + brii + c 







(5.8) 
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where 



An 2 , \2v 2 /p 2 4cr 2 /p 2 



— — Pi ~ vi , 6 = 2— — px - — u x , c = + pi — pi - 2—wi 



™2 VP2'^ "V ' ^2P2 V " 0"2 V ' ~ X ' r± \p2 r± P2~V 4n 2 ' 

(5.10) 

Now for our contour pi, a\ and v\ are always finite. For M\ and M 2 large and of the same 
order with |M 3 | << Mi,M 2 , the quantities p 2 and cr 2 given in (15.41) always remain large 
and positive, and the ratios (J 2 /p 2 , V2/P2, v 2 /a 2 etc. are bounded from above by finite 
numbers. Thus the quantities a, b and c remain finite. On the other hand in this limit 
we have 

p 2 a 2 -v%~ — X) (fMl + /3 2 M 2 2 ) + (A 2 + (1 - A) 2 + A(l - A)(a 2 + 5 2 )) M,M 2 . 

(5.11) 

In the range < A < 1 each term in this expression is non-negative, and (15.111) remains 
large and positive, - of order MiM 2 , - in the limit of large Mi, M 2 . This also shows that 
the matrix A defined in ( 15. 9p is nondegenerate in this limit and in fact has finite positive 
eigenvalues. As a result the term in the square bracket in (I5.8P is bounded from below by 
a finite number 

c-\{a b)A~ 1 (f \ , (5.12) 

and can never cancel the first term in (15.81) for any value of j and n\ for sufficiently large 
Mi and M 2 . Hence (15.81) cannot be satisfied. This shows that it is impossible to find a 
similtaneous solution to the eqs. (l5.1l) and (15.41) for n 2 7^ 0, and hence the contour C can 
be deformed to C" along the path (15.41) without encountering the poles of the integrand 
given in (15.11) for n 2 7^ 0. 

This establishes that in the limit of large charges the degeneracy remains the same 
in different domains in the moduli space up to nonperturbative terms. This result is 
consistent with the fact that for a black hole of charge (Q, P) in this theory we have 
a stable supersymmetric attractor for P 2 > 0, Q 2 > 0, P 2 Q 2 > (Q ■ P) 2 . Thus the 
near horizon geometry of these black holes is always given by this attractor point and is 
independent of the asymptotic moduli even if this requires the attractor flow to cross one 
or more walls of marginal stability. 
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6 Test of S-duality Invariance 



In §H we used S-duality invariance to determine the locations of the integration contour 
in the degeneracy formula in different domains in the moduli space. However this did not 
provide a test of S-duality. In this section we shall describe some tests of S-duality that 
one could perform. 

1. If there is an S-duality transformation that leaves the set Br invariant, then under 
such a transformation the contour C either should not transform, or should transform 
to another contour that is continuously deformable to C without passing through 
any poles. 

2. Analysis of [1] has shown that inside the left domain £ corresponding to the set of 
matrices Bl, the degeneracy is given by performing integration over the contour C 
described in ( I2.29p . Thus if there is an S-duality transformation that maps the set 
Br to the set Bl then such a transformation must map the contour C to the contour 
C or another contour deformable to C without passing through any pole. 

In fact for all values of N there is an S-duality transformation that maps Br to Bl. It is 
given by 



7 S J V° 1 

To see this we note that it maps the ( ^ ) e Br to ( \ .M 6 6i and ( \ ) E Bi 



^° ( n 1 ) G Bjj. Since two domains sharing two common boundaries must be identical 



l / 

the action of (I6.1I) must map Br to Bl- Using this transformation we can convert all tests 
of the first type into tests of the second type; all we need to do is to left multiply the 
duality transformation preserving Br by (I6.ip to construct a duality transformation that 
maps Br to Bl- 

We shall now try to verify that the transformation (I6.ip maps C to C or a contour 
deformable to C without passing through any poles. Using (14. 6ft . (14.191) we see that this 
transformation maps the contour C given in (12.151) to 

p 2 = M u a-. A 2 .!/, • - 2.Y.U :i . v 2 = NM 1 -M 3 , 

< pi < 1, < o x < N, < v x < 1 . (6.2) 
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This is different from C given in f!2.29j) for which (p 2 , o'2,v 2 ) = {Mi, M 2 , M 3 ). Thus we need 
to verify that we can deform the contour (16. 2h to C without encountering any pole. From 
the analysis in §(5]we already know that the poles at ( 15. II) for n 2 7^ are not encountered; 
thus we need to look for poles with n 2 = 0. They occur at 

(jv + n{5 — m{p + m 2 ) = , (6.3) 

with 

j 2 1 

nil E N 1Z>, rii, m 2 G j G 2 7L + 1, mini + — = - . (6-4) 
Taking the imaginary part of eq. (l6.3p we get 

jv 2 + n{o 2 - m{p 2 = . (6.5) 

For fixed j, n\, mi this describes a plane in the (p 2 , a 2 , v 2 ) space. Our job is to show that 
points in the (j) 2) '^2 ) v 2 ) space given in (12.291) and (16.21) lie on the same side of this plane 
so that we can deform them to each other without going through this plane. For this we 
need to show that the left-hand side of (I6.5p . evaluated at (12. 29ft and (16. 2ft have the same 
sign: 

(jM 3 + niM 2 - miMi) (j(NMi - M 3 ) + ni(N 2 Mi + M 2 - 2NM 3 ) - m x Mi) > . (6.6) 

We can simplify the analysis by setting Mi = M 2 = M. In this case the left hand side of 
(I6.6p has the form 

(rii — mi ){Nj + (N 2 + l)n a - mi}M 2 

+{(Nj + (N 2 + l)m - mi)] - (m - mi)(j + 2Nni)} M M 3 

-j(j + 2Nm)M 2 . (6.7) 

Using (16. 4p the coefficient of the M 2 term can be brought to the form 

fNj \ 2 9 9 N 2 , . 

[-± + m-mi\ +N 2 n\-—. (6.8) 

If n x ^ then this is strictly positive. If ni — then from (16 .4p we have j = ±1. We 
shall choose j = 1 by using the freedom of changing the signs of n x , mi, n 2 , m 2 and j 
simultaneously without changing the location of the pole. Since mi G N 'Z, in this case 
(16. 81) is strictly positive if mi 7^ 0, N, and vanishes for mi = 0, N. 
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Now as long as (16. 8p is non- vanishing and positive, we can make the first term of (16.71) 
dominate over others by taking M to be arbitrarily large, and hence (16. 7ft is positive as 
required. Thus we only need to worry about is the case n\ = 0, j = 1, mi = 0,N when 
the order M 2 term vanishes. In both cases (16.71) takes the form: 

N M M 3 — Mf . (6.9) 

Since M is large and positive, and M » M 3 > 0, this is strictly positive. This shows 
that (16. 6p holds for all j, mi, rti and hence we do not encounter any pole while deforming 
the contour to fl^29|) . 

For the special case of N — 1 one can consider another map that takes us from the 
set Br to El- It is 

a /3\ _ f 1 
7 Sj-{-l 

Under this map the vertices 0, 1 and ioo of the domain 7Z get mapped to the points ioo, 
— 1 and respectively. The latter set is precisely the vertices of the domain C. 

Using (14.61) . (14.191) we see that this transformation maps the original contour to [13] 

p 2 = M 2 , 3=2 = Mi, v 2 = M 3 , 
< pi < 1, < oi < N, < vi < 1 . (6.11) 

This exactly coincides with (12.291) up to an exchange of Mi and M 2 . Since exchange of Mi 
and M 2 does not change the integral (we could take Mi = M 2 ) we see that this constraint 
of S-duality is satisfied trivially. 

Returning to the case of general N, one can identify the following additional Ti(A^) 
transformation that maps the set Br to Bl'- 



(6.10) 



a (3\ _ ( 1 
7 5 ) ~ \-N 1 



(6.12) 



This maps the contour C to 



p 2 = Mi + M 2 - 2M 3 , a 2 = M 2 , v 2 = M 2 - M 3 , 

< pi < 1, < o x < N, 0<vi<l. (6.13) 

Consistency with duality invariance requires that we must be able to deform this contour 
to C without encountering any pole. One can proceed to analyze this exactly in the same 
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manner as we did for (16.21) and arrive at the condition that in order to be able to deform 
this contour to C given in (12. 29ft we need the following quantity to be positive (analog of 
eq.(EZj)): 

[rii — mi){ni — 2rrii + j)M 2 

+{j(ni - 2mi+j) + (m -mi)(2mi -])}MM 3 

+j(2m 1 -j)M 3 2 . (6.14) 

Using (16.4p the coefficient of the M 2 term can be brought to the form 

\ 2 1 

- + n 1 -m 1 j + m\ — -. (6.15) 

If mi 7^ it is strictly positive. For mi = from (16.41) we have j — 1. Since e in 
this case (I6.15P is strictly positive if ni ^ 0, —1, and vanishes for ni = 0, —1. As long as 
(16.151) is non-vanishing and positive, we can make the first term of (16.141) dominate over 
others by taking M to be arbitrarily large, and hence ( 16. 14[) is positive as required. Thus 
we only need to worry about is the case mi = 0, j = 1, rii = 0, —1 when the order M 2 
term vanishes. In both cases (16. 14H takes the form: 

MM 3 -M 2 . (6.16) 

Since M is large and positive, and M » M 3 > 0, this is strictly positive. This shows 
that (16.141) is positive for all j, mi, rii and hence we do not encounter any pole while 
deforming the contour (16 . 1 3[) to (12.291) . 
If we consider the Ti(N) element 

then it clearly leaves the set Br unchanged since it is given by a map from Br to Bl fol- 
lowed by the inverse of a map from Bl to Br. Our results establish that under this trans- 
formation the contour C is mapped to another contour that is continuously deformable 
to C. This is turn establishes that any power of go will also have the same property. 
We can now follow this up with the transformation (16.11) to construct a set of duality 
transformations that maps the set Br to Bl and maps the contour C to another contour 
deformable to C, thereby providing a test of the corresponding duality transformtion. For 
N = 1 this includes in particular the element ( 16.101) considered earlier. 
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Can every element of ^i(N) that preserves the set Br be expressed as a positive or 
negative power of g^l If so then our test of S-duality invariance of the degeneracy formula 
would be complete. To address this issue note that a Ti(iV) transformation that maps 
Br to Br must take adjacent walls to adjacent walls. Furthermore the map must be 
orientation preserving. Thus the action of any such transformation on Br must preserve 
the cyclic ordering of the walls and vertices. The map go indeed has this property. It 
moves the walls and vertices by two steps clockwise, taking to 1, ioo to 1 — i, 1 to 
1 — — etc. g^ 1 causes a shift by two steps in the anti-clockwise direction. Powers of 
go will move the walls clockwise or anti-clockwise by even number of steps. Are there 
elements of Fi(N) which shift the walls by odd number of steps? If so then by combining 
it with appropriate positive or negative powers of go we can generate a transformation 
that shifts every wall of Br by one step in the clockwise direction. Such a move will map 
to ioo, ioo to 1 and 1 to 1 — The unique SL(2,R) map that implements this is the 

VN -1/y/N' 



matrix ( ^ J . This is clearly not an element of Ti(N) for any A" other than 



N = 1. 

This would seem to indicate that all elements of Ti(N) preserving Br are generated 
by go. However there is an additional subtlety arising out of the fact^hat for A^ > 4 the 
element go has a pair of fixed points on the real line corresponding tc 



Too = -[l±Jl- — ] . (6.18) 




These provide accumulation points of the vertices of TZ; indeed if we begin with any 
vertex of 1Z and apply go or g^ 1 transformation successively to generate other vertices 
they accumulate at the points given in (16.181) . As a result there are infinite number of 
walls bordering 1Z, and a wall or a vertex that is situated within the range 




1 -\l 1 ~i-r I < r oo<^ I l + \/l-^7 I , (6-19) 



can never be related to a wall or a vertex lying outside this range by a go transformation. 
This opens up the possibility that there may be additional elements of ^i(N) which map 
the walls and vertices outside the range (16.191) to walls and vertices inside this range and 
vice versa, preserving the cyclic ordering. 



9 For N = 4 there is a single fixed point and hence it does not divide the vertices into two sets. Thus 
every pair of vertices separated by even number of steps can still be related by a go transformation. 
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To examine this issue we need to first identify some vertices lying within the range 
(16.191) . For this we shall first prove that the point 1/2 must be a vertex. If it is not a 
vertex then there must be a wall that goes over 1/2; by symmetry of the problem the 
vertices at the two ends of this wall must be situated symmetrically about 1/2. Let us 
take them to be 

l± P -- q -^. (6-20) 
2 q 2q 

with p,q relatively prime. This wall will correspond to a matrix \ ^ J with a/c = 



c d 

(q — 2p)/2q and b/d = (q + 2p)/2q. Now we see that if q is odd then (q ± 2p) are odd 
and hence both c and d must be divisible by 2. This of course is incompatible with the 
relation ad — be = 1. If on the other hand q is even (say 2m) then p must be odd and we 
can express (16.201) as 

(6.21) 

2m v ; 

Thus we have a/c = (m —p)/2m and b/d = (m + p)/2m. m must be odd so that (m ±p) 
are even; otherwise we again run into the problem of both c and d being even. Now since 
c and d cannot have a common factor, between m + p and m — p they should be able to 
cancel all the factors in 2m. This in particular will mean that their product (m 2 — p 2 ) 
must be divisible by m. Thus p 2 should be divisible by m. This is in contradiction with 
the assumption that p and q (and hence p and m) are relatively prime except for m = 1. 
The latter corresponds to the wall connecting and 1 and is not relevant for us. 

This shows that 1/2 must be a vertex (or accumulation point of vertices) of TZ. We can 
now begin with | and identify the vertices to the left and right of this using the analysis 
described below (13.331) . For odd N this gives three points 

i + (6.22) 

2 2N 2 2 2N K ' 

in ascending order. On the other hand for iV = 2M with M odd we get the three points 
to be 

11111 . 

, -, - + . 6.23 

2 2M' 2' 2 2M K ' 

If N is a multiple of 4 then one can show that there is no wall ending at 1/2 and hence 
1/2 must be an accumulation point. We shall not deal with this case here. We can now 
ask how go acts on these vertices. For N = 5 and N = 6 one can check that it moves 
the left-most vertex in the set (16.221) or (I6.23j) to the right-most vertex. By continuity we 
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can conclude that in the range (I6.19p . go moves every vertex and wall by two steps in the 

anti- clockwise direction. There is no contradiction with the fact that outside the range 

(16.191) go moves the points in the clockwise direction since the boundaries of this region 

are fixed points of go- Thus every wall in the range (I6.19P can be related by go action to 

one of the two walls connecting the vertices given in (16.221) or (16.231) for N = 5, 6. For 

N > 7 there are more than two vertices between a given vertex and its go image and the 

situation is more complicated. 

Restricting ourselves to the cases iV = 5, 6 we now ask the following question: is 

there a T\(N) transformation that maps the vertices lying inside the range (I6.19P to 

vertices outside the range ( 16.191) ? We focus on the possible action of this map on the 

three vertices given in (I6.22j) or (I6.23[) . If such a transformation exists then by left 

multiplying it with powers of go we can always bring the three final points to 1, zoo, 

or zoo, 0, 1/N. Examining the two cases separately we find that there is no such Ti(N) 

transformation mapping (I6.22|) or H6.23j) to these points (although for odd N there is a 

/ 2 -1 \ 

r (iV) transformation I , >, which takes the three points given in (16.221) 



N -(JV-l)/2 
to ioo, and 1/N). 

This finally establishes that for N < Q, all the elements of Ti(iV) which preserve the 
set Br are obtained by taking positive or negative powers of go- Since we have checked 
that these transformations take the contour C to another contour deformable to C without 
encountering any poles, this completes our test of S-duality invariance for N < 6. 
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